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Abstract 

The possibility to resolve narrow structures in reaction cross sections in calculations with the 
Lorentz integral transform (LIT) method is studied. To this end we consider a fictitious two- 
nucleon problem with a low-lying and narrow resonance in the 3 Pi nucleon-nucleon partial wave 
and calculate the corresponding "deuteron photoabsorption cross section" . In the LIT method the 
use of continuum wave functions is avoided and one works instead with a localized function ^f. In 
this case study it is investigated how far into the asymptotic region ^ has to be determined in 
order to obtain a precise resolution of the artificially introduced El resonance. Comparing with 
the results of a conventional calculation with explicit neutron-proton continuum wave functions 
it is shown that the LIT approach leads to an excellent reproduction of the cross section in the 
resonance region and of further finer cross section details at higher energies. To this end, however, 
for if one has to take into account two-nucleon distances up to at least 30 fm. 
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I. INTRODUCTION 

The LIT approach jl| allows the ab initio calculation of reaction cross sections, where 
a many-body continuum is involved. The great advantage of the LIT method lies in the 
fact that the knowledge of the generally complicated many-body continuum wave function 
is not required. In fact the scattering problem is reduced to a calculation of a localized 
function with an asymptotic boundary condition similar to a bound-state wave function. 
The LIT method has been applied to various electroweak cross sections in the nuclear mass 
range from A=3 to A=7. Among the applications are the first realistic ab initio calculations 
of the nuclear three- and four-body total photoabsorption cross sections 

aa 

as well as 

of the inelastic neutral current neutrino scattering off 4 He [4]. In addition first ab initio 
calculations were performed for the total photoabsorption cross sections of 4,6 He and 6 ' 7 Li 
with semirealistic forces {5], (3, B]. Other applications were carried out for the inelastic 
inclusive electron scattering cross section (see e.g. (si, Q). Further applications and a 
detailed description of the LIT method are presented in a recent review article 

In the past the LIT technique was mainly applied to cases where one has typically rather 
broad structures (quasielastic peak, giant dipole resonance). One may ask what is the 
power of the method in presence of narrow resonances in the continuum. In a recent LIT 
calculation of the (e,e') longitudinal and transverse form factors of 4 He 12| a resonance 
in the Coulomb monopole transition was observed but the width of the resonance could 
not be determined. Therefore in the present paper we want to address this problem and 
investigate how narrow resonances can be calculated with the LIT method. To this end we 
consider a fictitious two-nucleon problem with a low-lying and narrow resonance in the 3 Pi 
nucleon-nucleon partial wave and calculate the corresponding "deuteron photoabsorption 
cross section", which exhibits a pronounced El resonance at a photon energy of 2.65 MeV 
with a width of 270 KeV. 

The paper is organized as follows. In sect. II we give a short outline of the LIT calculation 
of the deuteron photodisintegration. In addition we discuss some aspects of the asymptotic 
behavior of the LIT solution \E'. In sect. Ill we describe the model for our fictitious NN 
interaction and show results for the corresponding 3 Pi phase shifts and the photoabsorption 
cross section calculated in the conventional way with a 3 Pi scattering wave function. The 
LIT results of our case study are discussed in sect. IV. 
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II. THE FORMALISM 



The deuteron total photoabsorption cross section is given by 

o*(u) = 4ir 2 aujR d (uj) , (1) 

where a is the fine structure constant, u is the energy of the photon absorped by the 
deuteron, and R d {ui) denotes the response function defined as 

R d (") = ^\(f\e\0)\ 2 S(co - E np - E d ) . (2) 

Here Ed and |0) are the deuteron bound state energy and wave function, E np and |/) denote 
relative kinetic energy and wave function of the outgoing np pair for a given two-nucleon 
Hamiltonian H, and O is the operator inducing the reaction (it is assumed that recoil effects 
are negligible). Here we consider the total photoabsorption cross section in unretarded dipole 
approximation, i.e. 

2 

e = 5>7f, (3) 

i=i 

where Zj and rf are the third components of position and isospin coordinates of the i-th 
nucleon. 

For a conventional evaluation of R d (u) one determines the np scattering wave functions 
for the induced NN partials waves ( 3 Po, 3 Pi, 3 P2 — 3 F2) and calculates the transitions matrix 
elements appearing in (JSJ). It is evident that for each of the three partial waves a separate 
cross section contribution can be defined. 

As already mentioned, with the LIT method one avoids the explicit calculation of scatter- 
ing wave functions. Here one proceeds in the following way. One first calculates the ground 
state wave function of the nucleus in question, in our case the deuteron. Then one has to 
solve the equation 

(H + E d -<r R -ia J )\*)=Q\0). (4) 

Since |\&) is localized one needs to apply only a bound-state technique for the solution of 
(jl]). It is convenient to perform multipole decompositions of left- and right-hand sides of 
(jlj). In the here considered deuteron case this leads to three separate equations, each one 
for a different (i=l,2,3). They correspond to the above mentioned three separate R d 
contributions due to the 3 Pq> 3 Pi, an d 3 -P2 — 3 F2 np final states. The calculation is carried 
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out for many values of an and a fixed 07. A theorem based on the closure property of the 
eigenstates of H shows that the sum of the overlaps (^l^j) of the three different multipole 
solutions of (jl]) corresponds to the Lorentz integral transform of R d , 

3 

L d (a R , 07) = X)<*l*> = / Rd ^) *R> a i) duJ > ( 5 ) 

i=l ^ 

where £ is a Lorentzian centered at ctr and with a width T = 207: 

C(u,a R ,a I ) = - \— — 2- (6) 

The parameters 0\r/j, for which (T4J) is solved, are chosen in relation to the physical problem. 
In fact <jj represents a kind of energy resolution for the response function, while the values 
of (Tr scan the region of interest. 

In a final step the transform (J5J) is inverted in order to obtain the response function 
and thus the cross section. In the following we use the standard LIT inversion method 
(alternative inversion methods can be found in 13(1)) however, with an extension that allows 
to take care of narrow resonances in the response function. 

The standard LIT inversion method consists in the following ansatz for the response 
function R d : 

M m 

ax 

R (-^np) ^ . ^■m~X_rn{,E n pi QJj) , (7) 



m=l 



where we have replaced the argument u of R by E np = uj — Ed. The Xm are given functions 
with nonlinear parameters aj. Here we take 

1/1 1 



Xi(E np , ai ) {Enp _ Ercs)2 + { i_ )2 ^ 1 + exp( _ 1) l + exp((K p -a 3 )/«3) 

Xm (E np , a t ) = E^ p exp(-^^) for m > 1 , (9) 

where in ((HI) we have set E res = oti and T = a 2 - It is evident that xi represents a resonance 
of Lorentzian shape, where the additional factor in brackets ensures that R d (E np ) is zero for 
E np = 0. For m > 1 the set \ m represents the usual basis function set for LIT inversions. 
Substituting such an expansion into the right hand side of (JSJ) one obtains 

ax 

L(0\R,07) = 22 C mXm{<?Ri<?IiUi) , (10) 
m=l 
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where 

~ / \ f jt-i Xm(E np ,Qii) 

Xm{<JR,(Tl,ai) = dE np — — — 2- ( n ) 

For given values of <X{ and M max the linear parameters c m are determined from a best fit 
of L d (<7R, 07) of (TTOT) to the calculated L d (0#,0/) of (J5j) for a fixed 07 and a number of 0# 
points much larger than M max . In addition one should vary the various nonlinear parameter 
a.i over a sufficiently large range. The parameter a^, however, can in general be determined 
from the known threshold behavior of the response function. In fact for our deuteron case 
one has transitions to P- waves of the final state and thus = 3/2. One starts the inversion 
by choosing a relatively low value of M max , e.g. M max = 6, then one selects the overall best 
fit and repeats the procedure for increasing M max up to the point that a stability of the 
inverted response is obtained and taken as inversion result. 

As pointed out above the key point of the LIT approach is the fact that one only works 
with a function, of finite norm. This localized function \l/ contains that information of the 
scattering process which is necessary to calculate reaction cross sections. The asymptotic 
behavior of ^ is described by an exponential fall-off with the argument ((2m/h 2 )ai) 1 ^ 2 , 
where m is the nucleon mass and 07 a parameter of the LIT method, which governs the 
resolution. The smaller 07 the better is the resolution and the more extended becomes 
In other words if one wants to resolve finer structures in the cross section one has to ensure 
that 07 is small enough and that ^ is calculated sufficiently far enough in the asymptotic 
region. 

In particular, in the following, we investigate up to which two-nucleon distance the lo- 
calized function \1/ has to be determined in order to obtain reliable results for the resonance 
cross section. In addition we also try to simulate the conditions of a LIT calculation for 
nuclei with more than two nucleons, where usually expansions on complete sets, e.g. hy- 
pherspherical harmonic functions, are used. In these cases the solutions \& are obtained in 
a kind of a spherical box of variable radius -R m ax- Beyond -R max ^ falls off rapidly. In order 
to simulate a similiar situation, in our case study we use an asymptotic boundary condition 
for ^ leading to a strong fall-off of ^ at a two-nucleon distance of R max and study the 
dependence of the results on i? max . 
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III. THE POTENTIAL MODEL 

In the following we will use the AV18 NN interaction [14J], but modify the 3 Pi potential in 
order to introduce a fictitious resonance. This is achieved by adding an attractive potential 
term, i.e. ^( 3 Pi) -> V( 3 Pi) + Kdd- We take 

57.6MeV , „,v / ,r - , . 

Kdd = (1 - exp(-2r 2 )) (l + exp(— )J (12) 

with the relative coordinate r in units of fm. 

In Fig. 1 we show the phase shift and in Fig. 2 the "deuteron photoabsorption cross 
section" cr^( 3 Pi) for the modified 3 Pi potential calculated with an explicit 3 Pi continuum 
wave function. The phase shift exhibits two resonances, one at E np =0A8 MeV and a second 
one at about 10.5 MeV. For <r^( 3 Pi) one finds a very pronounced and narrow resonance at 
cj=2.65 MeV with a width T of 270 KeV, while the second resonance is about four orders of 
magnitude weaker and has a width of about 5 MeV. There is a third cross section peak for 
cx^( 3 Pi) at 60 MeV, about 25 times stronger than the second resonance peak, with a rather 
large width, which cannot be ascribed to a 3 Pi resonance. 

IV. RESULTS AND DISCUSSION 

Before we investigate under which conditions one can reproduce the pronounced resonance 
of <t^( 3 Pi) of Fig. 2 via the LIT method, we first want to make another consideration. The 



smallest cross section width which has been resolved in previous L 

~ n 



T calculation is about 
Thus, in order to 



10 MeV (total photoabsorption cross sections of 6 He [6] and 2 H [ll|, 
better understand what was achieved in such calculations, let us first consider the deuteron 
case with the unmodified AV18 potential. Since we are particularly interested how far the 
asymptotic range of ^ has to be taken into account we proceed for the solution of (jlj) as 
follows. For a given value of P max we use a boundary condition that requires a strong fall-off 
of ^ at P max an d calculate the overlap (^|^)of ([5]) only in the range from r = to r = P max 
(for more details of such a calculation see [111 . Il6|). In Fig. 3 we show the conver gen ce of 
L d ( 3 Pi) with respect to P max taking, as for other previous LIT calculation of [Hi, Il5| . 
07 = 10 MeV. One sees that P max =10 fm is not sufficient, while P max =15 fm leads already 
to a rather good approximation of the final result. Convergence is essentially reached with 
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i? max =20 fm. One may ask whether 07=10 MeV is sufficiently small to resolve the peak 
structure of a d in the inversion. A minimal check can be performed in the following way. 
Let us assume that the cross section is a £-peak at E peak with size cr peak , i.e. 

a d (E np )=a^ k 5(E np -E pcak ). 

The resulting LIT is then given by the Lorentzian function 

L^(P pcak , a R , 7 ) = pcak £(P peak , a R , a f ) . 

In the right panel of Fig. 3 we show such an L d , taking as E p the peak position of 
L d (iCax = 30 fm) in the left panel of Fig. 3, in addition we set L^(P peak ) =L d (P peak ). It 
is readily seen that L d has a considerably larger width than L^ and thus should contain 
sufficient information of the peak structure for a reliable inversion. In fact in [ll] it was 
shown that a LIT calculation of L d with 0/ = 10 MeV leads precisely to the same cross 
section as a calculation with explicit np continuum wave functions. 

After these initial considerations we turn to a d ( 3 Pi) with the modified 3 Pi potential 
discussed in the previous section. We calculate L d ( 3 Pi) with various aj values, namely 
0/ = 10, 5, 2, 1, 0.5, and 0.1 MeV. First we discuss the convergence of the LIT with respect 
to -R max at low energies. In Figs. 4 and 5 we show the cases with 0/ = 1 and 0.5 MeV as 
examples. Different from the case considered in Fig. 3, P max =20 fm is not sufficient to reach 
a good convergence. For aj = 1 MeV one finds that the result for P max =30 fm differs up 
to about 1% from the converged LIT showing in addition a somewhat oscillatory behavior. 
With P max =40 fm one sees a much smoother curve and deviations are reduced to about 
0.5%. Very good convergence is reached with P max =50 fm with deviations of about 0.1%. 
A rather similar trend can be observed in Fig. 5 for 0/ = 0.5 MeV, but the deviations from 
the converged result are larger, e.g. for P max =30 fm one has deviations up to 8%, and the 
oscillatory behavior becomes more pronounced and is even visible for i? max =50 fm. For 
larger oj values (2, 5, 10 MeV) one finds a similar situation as shown in Figs. 4 and 5, but 
with increasing aj relative deviations become smaller and oscillations tend to vanish. For 
a j = 0.1 MeV it is not sufficient to take P m ax=80 fm in order to have a converged result. 
This is illustrated in Fig. 6. One observes a steady decrease of the peak height up to an 
-Rmax of almost 140 fm. In order to obtain a smooth result in the whole peak region one 
has to further increase P max up to 300 fm. It is interesting to observe in Fig. 6 how this 



convergence is reached with growing R mSLX . For values of 30, 60 and 90 fm two, three, and 
four separate Lorentzians, respectively, are visible. For even greater -R max the density of 
single Lorentzians increases further leading finally to a smooth curve for i? max =300 fm. The 
general picture looks very similar to LIT results obtained by solving \& via an expansion on 
a complete basis (see e.g. there an increase of basis states has a similarly effect on 

existing oscillations as a growing -R max in our case. 

As next point we want to check which value of 07 should be sufficient to obtain good 
inversion results at low energies, i.e. leads to a good description of the pronounced resonance 
peak in the cross section. For this purpose we perform a similar study as has been illustrated 
in the right panel of Fig. 3. In Fig. 7 we show the corresponding results for L d and L d for 
the various considered 07 values. For 07 = 10 MeV both results are essentially identical and 
for 07 = 5 MeV they are almost identical. A further reduction of 07 leads to more and more 
differences between L d and L d . Finally, for 07 = 0.1 MeV L d has a much larger width than 
L d . Thus one can be rather sure that a very good inversion result should be obtained for the 
pronounced low-energy resonance with our lowest 07 value, while 07=10 MeV should not be 
sufficient for a reliable inversion. For the other 07 values it is not clear beforehand whether 
they suffice, and in the following this will be investigated in detail. 

At first, in Fig. 8 we show the results of a minimal inversion, where in ([7j) we use M max =l 
and hence allow only for a single resonance in the cross section. It is evident that with 
decreasing 07 the resonance peak is steadily reduced and the resonance width steadily in- 
creased. In the right panel of Fig. 8 a comparison is made between the 07=0. 1 MeV result 
and the cross section calculated in the conventional way with a 3 Pi np continuum wave 
function. One finds an excellent agreement in the peak region. Thus it seems that it is 
necessary to use (i) large -R max values of far more than 100 fm and (ii) 07 values equal or 
smaller than the resonance width T. However, in the further discussion we will illustrate 
that, fortunately, this is actually not the case. 

Before we come to the discussion of more complete inversion results, in Fig. 9 we first 
show the LIT L d at higher energies considering three 07- and various -R max values. For 0=5 
MeV and -R m ax=30 fm one finds slight oscillations beyond 40 MeV, while a larger -R max leads 
to a smooth curve. The picture changes for 0=1 MeV, oscillations are present in all the 
three illustrated cases and start already beyond 20 MeV, it is also seen that the oscillation 
frequency (amplitude) increases (decreases) with growing i? max . For 0=0.1 MeV we only 
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show the result with i? max =300 fm. One notes that strong oscillations are present beyond 
20 MeV. 

As next point we turn to the full inversion of L d , i.e. not using as for the results of Fig. 8 
a restriction of M max to 1, but allowing for higher M max values. We first consider 07=1 MeV 
and take various R max values. With i? max =20 fm we do not obtain reasonable results, since 
the pronounced low-energy resonance cannot be reproduced, instead the best inversion fits 
are found for a vanishing width T of the resonance, which means that a 5-peak is preferred. 
As illustrated in Fig. 10 the situation changes for i? max =30 fm. For the low-energy region 
one has a very stable inversion result that starts already with M max =7 and leads to a very 
good description of the resonance peak. At higher energies the picture is different. None of 
the inversions describes the region of the second resonance correctly. At even higher energy 
inversions are improving with increasing M max but without reaching a very good agreement 
with the conventional calculation. The reason for not reproducing the high-energy cross 
section more precisely is due to the fact that the low-energy resonance contains almost all 
the El( 3 Pi) strength. Even for a high-energy a R the by far dominant LIT contribution 
stems from the low-energy resonance region, so that a precise inversion of the high-energy 
region requires also a very accurate calculation of the LIT. Similar results as for 07 = 1 
MeV are also found for 07 = 2 MeV, while our other 07 values do not lead to a sufficiently 
good description of the low-energy resonance if R max is taken equal to 30 fm. For -R max =50 
fm, also shown in Fig. 10, it does not come as a surprise that the low-energy region is 
desribed again in an excellent way. In addition for the higher M max values one finds a rather 
good description of the second resonance region and also the high-energy cross section is 
sufficiently accurate. In Fig. 11 we show results with -R max =80 fm and 07=0.5, 2, and 5 MeV. 
The smaller two 07 values lead to excellent descriptions of 0^( 3 Pi) in the whole considered 
energy range, only that with 07=0. 5 MeV the cross section is underestimated beyond 120 
MeV, while the 07=5 MeV result exhibits a to 5% too low resonance peak cross section. We 
do not show inversion results with 07=0.1 and 10 MeV. They are rather unstable and show 
oscillations with rather pronounced unphysical negative cross sections. 

We summarize our results as follows. We have performed a case study for a LIT calcu- 
lation with a pronounced resonance in the cross section. To this end the nucleon-nucleon 
interaction has been modified in the 3 Pi partial wave to generate a fictitious low-lying res- 
onance in the np continuum. The resulting "deuteron photoabsorption cross section" has 
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been calculated in two ways: (i) with an explicit 3 Pi np continuum wave function and (ii) via 
the LIT method, i.e. calculating the Lorentz integral transform with subsequent inversion 
in order to obtain the cross section. The conventional calculation leads to a cross section 
with three structures: (i) a very pronounced resonance with width of 270 KeV and a peak 
position of 0.48 MeV above the deuteron break-up threshold, (ii) a second resonance at 
10 MeV above threshold with a width of about 5 MeV, which is four orders of magnitude 
weaker than the first one, and (iii) a rather broad maximum at a photon energy of about 
60 MeV with a 25 times higher peak than the second resonance maximum. For a proper 
resolution of the dominant cross section structure, the pronounced low-energy resonance, 
one has to consider rather long-ranged solutions \& of the LIT equation, namely at least up 
to R max = 30 fm. The cross section results depend strongly also on the LIT parameter 07, 
which governs the resolution (the smaller 07 the better the resolution). It turns out that 
it is advantageous to work with 07 = 1 and 2 MeV, about 4-8 times larger than T, in this 
case -R max =30 fm is sufficient. Even larger 07 do not correctly reproduce the resonance 
cross section (07=5 MeV) or even lead to a 5-shape resonance (<T/=10 MeV). Very small 07 
values also lead to disadvantages, since a larger -R max is required, e.g. with 07=0.1 MeV, 
even taking R max = 300 fm the LIT is only poorly converged for energies beyond the first 
resonance region making inversions completely unreliable. Using the intermediate 07 values 
and setting i? max = 30 fm one obtains besides the excellent reproduction of the low-energy 
resonance also a rather good description of the broad maximum at higher energies, while 
the second, extremely weak, resonance is not well described. An increase of -R max to 50 
fm leads to a rather good description also of the tiny second resonance cross section, and 
with i? max = 80 fm there is further improvement so that an excellent agreement with the 
conventional calculation is achieved in the whole considered energy range if 07 is chosen 
between 0.5 to 2 MeV. 

One might think that for a good description of the low-energy resonance it could be 
sufficient to restrict the LIT inversion such that only a single resonance is allowed in the 
cross section. We have shown that this is actually not the case. Even taking -R max = 80 fm 
one obtains results which show an overestimation of the peak height and an underestimation 
of the width, e.g. with a/ = 5 MeV (1 MeV) one obtains a width of 90 KeV (180 KeV) 
instead of 270 KeV and a peak height which is overestimated by a factor of 2.7 (1.4). In 
addition the resonance is slightly shifted to higher energies. Only for 07 = 0.1 MeV and 
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-Rmax = 300 fm one finds a very good result for the resonance region. Fortunately, as pointed 
out in the previous paragraph, the picture changes if one performs a full LIT inversion, i.e. 
allows for structure besides the resonance. 

We think that these findings about the LIT resolution and the proper consideration of 
the asymptotic extension of the LIT solution ^ are generally important for LIT calculations. 
In particular for methods where ^ is expanded on a complete set it has to be ensured that 
the solution extends far enough into the asymptotic region if one wants to properly resolve 
cross section structures with a small width. If this is guaranteed the LIT approach opens 
the possibility to calculate also narrow cross section resonances with great precision. 
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FIG. 2: Photodisintegration cross section cr^( 3 Pi) with modified V( 3 P\) potential at low energy 
(left) and for extended energy range (right). 
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FIG. 3: Lorentz integral transform L d ( 3 Pi) with AV18 NN potential in low-energy region for 
a i = 10 MeV and various values of R m & x (left) and comparison of the R m & x = 30 fm result with 
Lj (right). 
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FIG. 4: Lorentz integral transform L d ( 3 Pi) with modified V( 3 P\) potential in first resonance 
region for aj = 1 MeV and various values of i? max (top) and ratio R=L d (i? max )/L d (i? max = 80 fm) 



(bottom). 
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FIG. 5: As Fig. 4 but for 07 = 0.5 MeV 




FIG. 6: As Fig. 4 but for 07 = 0.1 MeV and various values of R max in both panels. 
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FIG. 7: Lorentz integral transform L d ( 3 Pi) (full curves) with modified V( 3 Pi) potential in first 
resonance region in comparison to L d (dashed curves) for various values of 07 with i? max = 80 fm, 
except for 07 = 0.1 MeV where R mSLX is equal to 300 fm. 
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FIG. 8: Cross section a^( 3 P\) with modified V( 3 P\) potential in first resonance region obtained 
from "minimal" inversion with M max = 1 (see text) for various 07 values (left) and comparison 
of the 07 = 0.1 MeV result (full curve) with result of conventional calculation with explicit np 
continuum wave functions (crosses). 
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FIG. 9: Lorentz integral transform L d ( 3 Pi) with modified V( 3 Pi) potential beyond first resonance 
region for 07 = 5 MeV (top), 1 MeV (middle), and 0.1 MeV (bottom) and various values of i? m ax- 



17 



12 
10 

! 8 

V 4 
2 

10 
10 

! 8 

V 4 
2 




1 


i 1 i 

< M 

' max 

\ ■■■ 8 


/ 


\ --- 15 


- / 


\ — 30 




\ R =30 fm 

\ max 




\ o=l MeV 


- / 




1 *- 


— 


i M 

X max 

A ■ ■ ■ 7 


/ 


\ --- 15 




\ — 30 




\ R =50 fm 

\ max 




\ o=lMeV 
1 ^ *- 





10"' 



10" 



10 




"max 15 

o=lMeV — 30 




3.5 



to [MeV] 



20 40 60 80 100 120 140 
co [MeV] 



FIG. 10: Cross section cj^( 3 Pi) with modified V( 3 P\) potential obtained from inversion of 
L d ( 3 Pi,cr/ = 1 MeV, i?max) with various values of M max and for i? max = 30 fm (top) and 50 
fm (bottom) in first resonance region (left) and at higher energies (right); also shown results of 
conventional calculation with explicit np continuum wave functions (crosses). 
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FIG. 11: As Fig. 10 but for i? max = 80 fm and a T = 5, 2, and 0.5 MeV. 
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